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Abstract. In this talk I review studies of hadron properties in bosonized chiral quark models for the 
quark flavor dynamics. Mesons are constructed from Bethe-Salpeter equations and baryons emerge 
as chiral solitons. Such models require regularization and I show that the two-fold Pauli-Villars 
regularization scheme not only fully regularizes the effective action but also leads the scaling laws 
for structure functions. For the nucleon structure functions the present approach serves to determine 
the regularization prescription for structure functions whose leading moments are not given by matrix 
elements of local operators. Some numerical results are presented for the spin structure functions. 



1. Introduction 

In this talk I review investigations of hadron properties in the Nambu-Jona-Lasino (NJL) 
model [1]. This is a particularly simple model for the quark flavor interactions with the 
important feature that the quarks can be integrated out in favor of meson fields [2]. The 
resulting effective action for these mesons possesses soliton solutions [3]. According to 
the large-A^c picture [4] of Quantum-Chromo-Dynamics (QCD) these solutions are in- 
terpreted as baryons. 

The construction of hadron wave-functions is not possible in QCD. This represents a 
main obstacle for the computation of hadron properties from first principles. As the NJL 
model adopts the symmetry properties of QCD, the current operators in the model cor- 
respond to those of QCD. As a consequence, matrix elements of the current operators as 
computed in the model are sensible and their comparison with experimental data is mean- 
ingful. In particular, it is interesting to analyze the hadronic tensor that parameterizes the 
deep-inelastic-scattering (DIS) and confront the model predictions with empirical data. 
This picture has led to interesting studies of hadron structure functions in bosonized chi- 
ral quark models. Here I will present the results of refs. [5-7]. These studies build up a 
consistent approach by computing the hadronic tensor (or equivalently the forward virtual 
Compton amplitude) from the gauged meson action. For the nucleon structure functions 
similar studies have been reported in refs. [8-10]. There no attempt to compute the struc- 
ture functions from the gauged action was made but rather it was assumed that the model 
predictions for the constituent quark distributions can be identified with QCD quark distri- 
butions. I refer to those articles for a more expatiated presentation of numerical results. In 
addition, I refer to the review articles [3] for comprehensive discussions of model predic- 
tions for static baryon properties such as magnetic moments, axial charges or the hyperon 
spectrum. 
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This talk is organized as follows. In Section 2 I introduce the NJL model as an effective 
meson theory and utiUze pion properties to determine the model parameters. Section 3 de- 
scribes the subtleties for extracting the structure functions that arise in this model from reg- 
ularization. The pion structure function is considered as an example. In Section 4 I review 
the construction of baryon states in the soliton picture. The following Section sketches 
the computation of nucleon matrix elements of the hadronic tensor and the extraction of 
the structure functions in the Bjorken limit. Finally in Section 6 I present some numerical 
results for the spin structure functions gi and 52 and compare them to experimental data 
by means of the transformation to the infinite momentum frame and subsequent DGLAP 
evolution. Section 7 serves as a short summary. 



2. The NJL Model for Chiral Dynamics 

The NJL model is a quark model with a chirally invariant quartic quark interaction. 
Bosonization is achieved semiclassically by introducing effective meson fields for the 
quark bilinears in that interaction. Then the quark fields are integrated out by functional 
methods. This yields an effective action for meson degrees of freedom, 

A[S, P] = -iNcTrj^log [t^ -{S + ij^P)] -^Jd^x tvViS, P) . (1) 

Here V is a local potential for the effective scalar and pseudoscalar fields S and P, 
respectively, that are matrices in flavor space. In the NJL model the potential reads 
V = 5^ + + 2rhoS with mo being the current quark mass matrix. Since the in- 
teraction is mediated by flavor degrees of freedom, the number of colors, Nc, is merely a 
multiplicative quantity. The functional trace (Tr) originates from integrating out the quarks 
and induces a non-local interaction for S and P. For simplicity I will only consider the 
isospin limit for up (u) and down (d) quarks: mo,u = mo,d = mo. 

A major concern in regularizing the functional (1), as indicated by the cut-off A, is to 
maintain the chiral anomaly. This is achieved by splitting this functional into 75-even and 
odd pieces and only regulate the former. 
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iT> = i^-{S + i^nP) and iDs = -i^ - {S - i-y^P) . 
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The double PauU-Villars regularization renders the functional (1) finite with co = 1, Ao 



0, e:=o 



. The 75 -odd piece is conditionally finite and not regularizing it, repro- 



duces the chiral anomaly properly. For sufficiently large G one obtains the VEV, {S) = m1l 
that parameterizes the dynamical chiral symmetry breaking from the gap-equation. 
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Substituting S = (S) = m1l in eq. (1) shows that m plays the role of a mass and is 
therefore called the constituent quark mass. 

In the next step I utilize pion properties to fix the model parameters and introduce the 
isovector pion field tt via 

S + iPjs = m (Uy^ = m exp f i — 75 n ■ t) . (5) 

\ m / 

Sandwiching the axial current between the vacuum and a single pion state yields the pion 
decay constant /^r = 93MeV in terms of the polarization function Il{q'^,x), 

fn = 4:Ncmg / dxll{ml,x) 
Jo 

x) = ^ c„ — — [k^ + x{l - x)ml - - + ie] , (6) 

where m^r = 138MeV is the pion mass. The Yukawa coupling constant, g, is determined 
by the requirement that the pion propagator has unit residuum, 

? = [ '^'^ Kn(rri^, x)] . (7) 

In the chiral limit (rn,^ = 0) this simplifies to = m/g. Finally the current quark mass is 
fixed from the condition that the pole of the pion propagator is exactly at the pion mass, 

mo=^NcmGm\ [ dxn{ml,x). (8) 
Jo 

It is also worthwhile to mention that expanding eqs. (2) and (5) to Unear and quadratic order 
in TT and u^, respectively, yields the correct width for the anomalous decay — > 77. This 
is the direct consequence of not regularizing the 75 -odd piece. 

Before discussing nucleons as soUtons of the bosonized action (1) and the respective 
structure functions it will be illuminating to first consider DIS off pions. 



3. The Compton Tensor and Pion Structure Function 

DIS off hadrons is parameterized by the hadronic tensor W^'^ {p, q) where q is the momen- 
tum transmitted from the photon to the hadron with momentum p. 

The tensor W^" (p, q) is obtained from the hadron matrix element of the current commu- 
tator by Fourier transformation and is parameterized in terms of form factors that multiply 
the allowed Lorentz structures. These form factors are obtained by pertinent projection 
of the hadronic tensor. Finally the structure functions are the leading twist contributions 
of the form factors. These contributions are obtained from computing W^" (p, q) in the 
Bjorken limit: = — — > 00 with x = Q^/p ■ q fixed. That is, subleading contributions 
in l/Q^ are omitted. 

An essential feature of bosonized quark models is that the derivative term in (1) is for- 
mally identical to that of a non-interacting (or asymptotically free) quark model. Hence 



the current operator is given as = qQ'y'^q, with Q a flavor matrix. Expectation values 
of currents are computed by introducing pertinent sources in eq. (2) 

iT> — > iD + and iT>5 — > iDs - (9) 

and differentiating the gauged action (1) with respect to v^. In bosonized quark models it 
is convenient to start from the absorptive part of the forward virtual Compton ampUtude^ 

T'^'^ip, q) = J rf^Ce'«-« (P, s\T {J'^iOrm \p, s) , W'^'^ip, q) = [T'^'^ip, q)] 



because the time-ordered product is straightforwardly obtained from 



(10) 



(11) 



as defined from eq. (2) with the substitution (9). 

Pion-DIS is characterized by a single structure function, F{x). For its computation 
the pion matrix element in the Compton amplitude (10) must be specified. For virtual 
pion-photon scattering it is obtained by expanding eqs. (2) and (5) to second order in 
both, TT and Vfi. Due to the separation into D and D5 this calculation differs considerably 
from the simple evaluation of the 'handbag' diagram. For example, isospin violating and 
dimension-five operators appear for the action (2). Fortunately all isospin violating pieces 
cancel yielding 

F{x) = •i{4Nc9')^ [mlU{ml,x)] , 0<x<l. (12) 

The same result is obtained by formal treatment of the divergent handbag diagram and ad 
hoc regularization [1 1]. The cancellation of the isospin violating pieces is a feature of the 
Bjorken limit: insertions of the pion field on the propagator carrying the infinitely large 
photon momentum can be safely ignored. Furthermore this propagator can be taken to 
be the one for non-interacting massless fermions. This impUes that also the Pauli-Villars 
cut-offs can be omitted for this propagator. That, in turn, leads to the desired scaling 
behavior of the structure function in this model and is a particular feature of the Pauli- 
Villars regularization. A priori it is not obvious for an arbitrary regularization scheme that 
terms of the form Q^/A^ drop out in the Bjorken Umit. 

From eqs. (7) and (12) it is obvious that F{x) = 5/9 for < .t < 1 in the chiral limit 
(TOtt = 0). It must be noted that this refers to the structure function at the (low) energy 
scale of the model. To compare with empirical data, that are at a higher energy scale, 
the DGLAG program of perturbative QCD has to be applied to F{x) to include the InQ^ 
corrections. Such studies [12] show good agreement with the experimental data for F{x). 



^The momentum of the hadron is called p and its spin eventually s. 



4. The Nucleoli as a Chiral Soliton 



Solitons are a non-perturbative stationary configurations of the meson fields. To determine 
that configuration for the meson theory (1)1 consider the hedgehog ansatz 

Uuir) = exp {if ■ fF{r)) and (Uii{f)y^ = exp {ij^f ■ fF{r)) 

for the pion field (5). The corresponding single particle Dirac Hamiltonian reads 

h = a ■ p + (3m [cosF + i^^f ■ f sinF] . (13) 

Evaluating the action functional (2) in the eigenbasis of h gives the energy functional in 
terms of the eigenvalues, Cq, [13] 

E[F] = ^{1- sign(ev)) - ^ E E ^4 - \/ + k] 

+mlfl j d\ (1 - cosF) (14) 

for a baryon number one configuration. Here V denotes the unique quark level that is 
strongly bound by the soliton. Its explicit occupation takes care of the total fermion number 
and thus this level is referred to as the valence quark. It should not be confused with the 
valence quarks in the parton model. Furthermore eL°^ are the eigenvalues of /i^^^ = a ■ p + 
P m . The soliton profile F{r) is then obtained from extremizing E self-consistently [3]. 

States possessing good spin and isospin quantum numbers are generated by taking the 
zero-modes to be time dependent [14] 

U{f,t) = A{t)UYi{r)A\t) , (15) 

which introduces the collective coordinates A{t) G SU{2). The action functional is ex- 
panded [15] up to quadratic order in the angular velocities 

if ■^^2A\t)A{t) . (16) 

The coefficient of the quadratic^ term defines the moment of inertia^, Q:^[i^]. Nucleon 
states \N) are obtained by canonical quantization of the collective coordinates, A{t). This 
is analogous to quantizing a rigid rotator and allows to compute matrix elements of opera- 
tors in the space of the collective coordinates [14]: 

{N\\tv{Tc,A^nA)\N) = -l{N\hJb\N) and Vt = -J/a''[F\, (17) 
where and denote isospin and spin, respectively. 



liner term does not arise due to isospin symmetry. 
^Functional integrals are evaluated using the eigenfunctions 0a of the Dirac Hamiltonian (13) in 
the background of the chiral angle F{r). Thus all quantities - like the moment of inertia - turn into 
functionals of F{r). 



For later use I note that the valence quark wave-function receives a first order cranking 
correction 

*v(f ,i) = e-^-*A(0 I Mr) + I E 0,(r) ^^'^-^l^^ I , (18) 

where (jy^ (f ) are the eigenfunctions of /i in eq. (13). The moment of inertia, [F] is order 
Nc, thus, upon quantization (17), this rotational correction is subleading in l/Nc- 



5. Nucleon Structure Functions 

DIS off nucleons is described by four structure functions: Fi (x) and F2 (x) are insensitive 
to the nucleon spin while the polarized structure functions, gi (x) and g2{x), are associated 
with the components of the hadronic tensor that contain the nucleon spin. 

As argued in section 3, the quark propagator with the infinite photon momentum should 
be taken to be the one for free and massless fermions. Thus, it is sufficient to differentiate 
(Here D and D5 are those of eq (3), i.e. with Vfj, = 0.) 



with respect to the photon field V/^. I have introduced the (. . .)5 description 



(19) 



Sn 



to account for the unconventional appearance of axial sources in D5, cf. ref. [7]. Substitut- 
ing eq. (15) for the meson fields that are contained in D and D5, computing the functional 
trace up to subleading order in 1/A''c using a basis of quark states obtained from the Dirac 
Hamiltonian (13), yields analytical results for the structure functions. 1 refer to ref. [7 J for 
detailed formulas for other structure functions and the verification of the sum rules that 
relate integrals over the structure functions to static nucleon properties. As an example I 
restrain myself to list the contribution to gi {x) which is leading order in 1 /Nc'- 
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(20) 



where the subscript (P) indicates the pole term. 

Before discussing numerical results I would Uke to mention the unexpected result that 
the structure function entering the Gottfried sum rule is related to the 75-odd piece of the 



action and hence does not undergo regularization. This is surprising because in the parton 
model this structure function differs from the one associated with the Adler sum rule only 
by the sign of the anti-quark distribution. The latter structure function, however, gets 
regularized in the present model, in agreement with the quantization rule for the collective 
coordinates that correspond to the isospin operator that involves the regularized moment 
of inertia, o? . 

6. Numerical Results for Nucleon Structure Functions 

Unfortunately numerical results for the full structure functions in the double Pauli-Villars 
regularization scheme, i.e. including the properly regularized vacuum piece are not yet 
available. However, in the Pauli-Villars regularization the axial charges are saturated to 
95 % or even more by their valence quark (18) contributions once the self-consistent soliton 
is substituted. This provides sufficient justification to consider the valence quark contribu- 
tion to the polarized structure functions as a reliable approximation since e.g. the zeroth 
moment of the leading structure function g\ is nothing but the axial current matrix element. 
This valence quark level is that of the chiral soliton model and, as already mentioned, its 
contributions to the structure functions should not be confused with valence quark distribu- 
tions in parton models. In general, it should be stressed that the present model calculation 
yields structure functions, i.e. quantities that parameterize the hadronic tensor, but not 
(anti)-quark distributions. The latter would require the identification of model degrees of 
freedom with those in QCD. However, here only the symmetries (namely the chiral sym- 
metry) and thus the current operators in the hadronic tensor are identified. 

As in the case for the pion, the model calculation yields the nucleon structure function 
at a low energy scale. In addition the soliton is a localized object. Thus the computed 
structure functions are frame-dependent and one frame has to be picked. The appropri- 
ate choice is the infinite momentum frame (IMF) not only because it makes contact with 
the parton model but also because it is that frame in which the support of the structure 
functions is limited to the physical regime < a; < 1. Choosing the IMF amounts to the 
transformation [16,17] 

/iMF(a;) = ^ /RF(-ln(l - x)) , (21) 

where fRF{x) denotes the structure function as computed in the nucleon rest frame. So the 
program is two-stage, first the transformation of the model structure function to the IMF 
according to eq (21) and subsequently the DGLAP evolution program [18] to incorporate 
the resumed InQ^ corrections. In the current context it is appropriate to restrain oneself 
to the leading order (in Ug) in the evolution program because higher orders require the 
identification of quark and antiquark distributions in the parton models sense. In the present 
model calculation this is not possible without further assumptions'*. The low energy scale, 
Qq = 0.4GeV^, at which the model is assumed to approximate QCD has been estimated in 
ref. [5] from a best fit to the experimental data of the unpolarized structure function, F2 {x). 



''We assume, however, that the gluon distribution is zero at the model scale. 
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Figure 1. Model predictions for tlie polarized proton structure functions xgi 
(left panel) and xg2 (right panel). The curves labeled 'RF' denote the results as 
obtained from the valence quark contribution to (19). These undergo a projec- 
tion to the Infinite momentum frame 'IMF' (21) and a leading order 'LO' DGLAP 
evolution [18]. Data are from SLAC-E143 [21]. 

The same boundary value is taken to evolve the model prediction for polarized structure 
function, gi (x), in the IMF to the scale of several GeV^ at which the experimental data 
are available. For the structure function g2{x) the situation is a bit more compUcated. First 
the twist-2 piece must be separated according to [19] 
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(x) = -gi{x) + [ — 

Jx y 



9i{y) 



(22) 



and evolved analogously to gx {x) (which also is twist-2). The remainder, 52 {x) —gY^ {x), 
is twist-3 and is evolved according to the large- A^c scheme of ref. [20]. Finally, the two 
pieces are again put together at the end-point of the evolution, Q^. In figure 1 1 compare the 
model predictions for the linearly independent polarized structure functions of the proton 
to experimental data [21]. In figure 2 1 compare the model predictions for both the proton 
and the neutron (in form of the deuteron) not only to the recently accumulated data but also 
to other model predictions. Surprisingly the twist-2 truncation, i.e. eq. (22) with the data 
for gi (x) at the right hand side, gives the most accurate description of the data. However, 
also the chiral sohton model predictions reproduce the data well. Bag model predictions 
have a less pronounced structure. 

Recently, precise data [25] have become available for the neutron asynnmetry 



Fi{x,Q^) 



(23) 



It is therefore challenging to study this quantity in the present model. As subleading twist 
contributions are omitted, this amounts to computing the ratio gi{x,Q'^)/Fi{x,Q'^), for 
the neutron. The resulting ratio is shown in Fig. 3 together with data. It is interesting to 
note that while the ratio at the model scale, Qo, becomes large and negative at small x, the 
DGLAP evolution causes it to bend around so that it actually tends to zero as a; ^ 0. This 
behavior is also observed from the data, as is the change in sign at moderate x. The position 
(x « 0.25) at which this change occurs seems somewhat lower than the preliminary JLab- 
data [25] suggest and insensitive to the end point of evolution. Once evolution has set in at 
a moderate point Q^, the evolution to even higher has insignificant effect. 



0.04 
0.02 



O.OO 

— 0.03 

— 0.04 

X — 0.06 
0.04 

0.02 

O.OO 
— 0.03 
— 0.04 

0.02 0.05 O.IO 0.20 0.50 l.OO 

X 

Figure 2. Model predictions for the polarized proton structure functions xg2 for 
proton and neutron (deuteron) and comparison with data from E143 [21] (open 

diamond) and E155 [22] (open square) and their combination (solid circle). The 
full line is the twist-2 truncation (22) of data for gi{x). Dashed-dotted [23] and 
dotted [24] lines are bag model calculations, the short dashed lines represent 
the present chiral soliton model [6] and long dashed line that of ref. [10]. (This 
is a slightly modified figure from ref. [22].) 

7. Conclusions 

I have discussed a chiral quark model for hadron phenomenology. In particular, I consid- 
ered the bosonized NJL model as a simplified model for the quark flavor dynamics. Al- 
though the bosonized version is a meson theory, the quark degrees of freedom can indeed 
be traced. This is very helpful for considering structure functions. Additional correlations 
are introduced due to the unavoidable regularization which is imposed in a way to respect 
the chiral anomaly. Hence a consistent extraction of the nucleon structure functions from 
the Compton amphtude in the Bjorken limit leads to expressions that are quite different 
from those obtained by an ad hoc regularization of quark distributions in the same model. 
I also showed that within a reliable approximation the numerical results for the spin depen- 
dent structure functions agree reasonably well with the empirical data. 
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